function fa=FA(tc, tv, h)
    fa = abs(1 + h*(tc - tv));
